Exponential Temperature Dependence of Penetration Depth in MgB2 
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The temperature dependence of the London penetration depth, A(r), was measured in polycrys- 
taUine MgB2 samples by a high-resolution, radio frequency technique. A clear exponential temper- 
ature dependence of A(T) was observed at low temperature, indicating s-wave pairing. A fit to the 
data gives an effective energy gap A of 2.8 ± 0.4 meV (2A/fcsTc = 1.7 ± 0.2), which is significantly 
smaller than the standard BCS weak coupling value of 3.5. This may be explained either by an 
anisotropic gap or the influence of low frequency phonon modes. 
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The recent discovery Q of superconductivity at 39 K 
in the binary compound MgB2 has sparked a large num- 
ber of investigations into its physical properties. A cru- 
cial question is whether its high Tc can be explained by 
a phonon mediated pairing interaction within the usual 
BCS-Eliashberg framework. A first step in answering this 
question is to determine the symmetry of the supercon- 
ducting order parameter and the nature of the low energy 
excitations. 

Several experiments so far, have given strong indica- 
tions that the superconductivity is of the s-wave BCS 
type. An important experiment has been the observa- 
tion of a large boron isotope effect In conventional 
superconductors, tunneling measurements played a key 
role in verifying the BCS theory, however, in the cuprate 
superconductors this probe has not proved quite as use- 
ful because the short coherence length ^ makes junction 
fabrication difficult. The majority of tunneling measure- 
ments on MgB2 so far show a s-wave like gap but the 
values of the deduced gap vary widely between 2 and 7 
meV 1^-^ . The magnetic penetration depth A has proven 
to be a powerful probe of the gap symmetry in both classi- 
cal and high Tc superconductors and gives unambiguous 
information about the nature of the low energy excita- 
tions in superconductors. As the length scale probed is of 
order A rather than ^, for high k materials, penetration 
depth measurements are much less sensitive to surface 
quality than tunneling, and hence we expect the results 
to be more representative of the bulk. 

There have been several recent reports of measure- 
ments of A(T) for MgB2 by ac-susceptibility muon 
spin rotation 0| and optical conductivity techniques. 
These authors conclude that at high temperature A(T) 
follows a power law dependence [A(T) ^ (Refs. [0,D) 
and A(r) - T^ "^ (Ref. g)], which seems to be at odds 
with tunneling and other measurements which indicate 
that there is a sizeable s-wave gap. In this paper we 
present high resolution measurements of the temperature 
dependence of the penetration depth in polycrystalline 
samples of MgB2. We find strong evidence for a predom- 
inately exponential temperature dependence of A at low 
temperature consistent with s-wave BCS behavior. The 



gap deduced from fits to the data however, is significantly 
smaller than the BCS weak coupling value. This may be 
explained by either a strong anisotropy of the gap or the 
presence of low energy phonon modes. 

Three different types of polycrystalline MgB2 samples 
were studied. Samples A and C, were produced from 
commercially available (Alfa Aesar) powder. Sample A 
was produced by casting the powder directly in epoxy 
resin (weight ratio MgB2:Epoxy 1:6.5). The commer- 
cial powder has a relatively wide range of grain sizes 
which complicates quantitative analysis of the data. For 
the analysis to give a reliable estimate of the absolute val- 
ues of A(T) grain sizes ^ 5A are required ||l^. To obtain 
such a distribution the commercial powder was ground 
in an agate mortar and then sedimented in acetone for 1 
hour |l^. The resulting powder was then cast in epoxy 
(~6% by volume) to produce sample C. The grain size 
distribution was measured from scanning electron micro- 
scope images. The diameter of 96% of the grains was less 
than 1/im, and the mean diameter was 0.56/im. Sam- 
ple B was fabricated by reacting boron powder and Mg 
flakes in an Argon atmosphere at high temperature. This 
procedure produced a dense polycrystalline pellet. X- 
ray diffraction and scanning electron microscopy showed 
that as grown samples of type B contained some resid- 
ual flakes of unreacted Mg. Initial measurements found 
that this changed the low temperature behavior of A(r) 
markedly |jl^. To remove this contamination we im- 
mersed these samples in a dilute solution (~ 0.5 %) of 
HCL and ethanol. Tc was unaffected by this. Previous 
studies have shown that the commercial MgB2 powder 
is 98% phase pure, the remaining 2% being mostly com- 
posed of iron, carbon and oxygen |l^ . The susceptibility 
of these impurities is too weak to effect the results pre- 
sented here [|l|,|l|l 

Measurements of penetration depth were performed in 
a tunnel diode oscillator operating at 11.9 MHz with 
frequency stability of a few parts in 10^° Hz" 2. For the 
powder samples in the present work this translates to a 
resolution in A of 10^^'^m. A particular feature of our ap- 
paratus is the very low value of the ac-probe field which 
we estimate to be ~ l/iT. Ambient dc fields are shielded 
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to a similar level with a mu-metal can. Changes in the 
oscillator frequency are directly proportional to the in- 
ductance of the probe coil and hence to the susceptibility 
of the sample. For single crystal samples, this frequency 
shift may be directly related to changes in the penetra- 
tion depth using the known sample dimensions fl^ ]. In 
general for polycrystalline samples the relation between 
the measured susceptibility (per unit volume of super- 
conductor) X ^iid A is more complicated and depends on 
the size distribution of the grains. For well separated 
grains 



E(i-f coth(5^) + 



3A^ 



(1) 



where Ni is the measured number of grains of radius 
[ pot . For sample C, A(T) was determined from the mea- 
sured x{T) and grain size distribution by solving Eq. (|l|) 
at each point. For samples A and B the wide distribu- 
tion of grain sizes and (for sample B) the possibility of 
interconnection between the grains makes a quantitative 
determination of A(T) impossible. We note however, that 
if A(T) becomes temperature independent, then so will 
x{T), irrespective of the detailed dependence of x{T) on 
A and r^. 

In Fig. we show the temperature dependence of x for 
the three different samples. For samples A and B, x was 
normalized to —1 at low temperature, whereas for sample 
C, X calculated from the measured volume of super- 
conductor and is thus in absolute units. Ax(T) denotes 
the change in x relative to its value at our base temper- 
ature {T—IA K). It is apparent that in all three samples 
Ax{T) becomes essentially temperature independent be- 
low around 5 K. This is the main result of this paper 
and is a clear indication of an finite energy gap in all 
directions in fc-space. A large distribution of grain sizes 
or anisotropy of A will effect the detailed T dependence 
of X (especially at higher temperature) but cannot cause 
such an abrupt flattening off of Ax(T). 

The material parameters of MgB2 deduced from other 
measurements put it firmly in the clean, local limit 
(A(0) ~ 0.2/im , ^(0) -50 A and mean free path £ - 600 
A p7| , p^ ), so Ax{T) is directly proportional to changes 
in the London penetration depth (Al). For a s-wave su- 
perconductor, at sufhciently low temperature, irrespec- 
tive of the coupling strength, Xl{T) should follow the 
BCS behavior 



AA(r) ~ A^ 



ttAq 
2kBT 



exp 



-Aq 
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were Aq is the value of the energy gap at zero tempera- 
ture and Ag = A(0) in the weak coupling case. In Fig. |l| 
we show a fit of Eq. (||) to the data for the three sam- 
ples. The fits give consistent values for Aq which average 
33 ± 4 K (the uncertainty reflects the different values of 
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FIG. 1. Temperature dependence of the susceptibility (x) 
of samples A, B and C. At low temperature (main panels) Ax 
is directly proportional to the temperature dependent pene- 
tration depth AA(T). The solid line is a fit to the s-wave 
BCS form [Eq. (^)], and the dashed line to a behavior. 
The insets show x(2^) over the full temperature range. 



Aq between the three samples and as the upper limit of 
the fit (T^ix) was varied from 15 K to 20 K). It should 
be noted that Eq. (|^) is not a perfect fit to the data and 
particularly at very low temperature where the measured 
AA(T) has a larger temperature dependence that would 
be indicated by Eq. (||). If taken literally this may indi- 
cate the presence of a small amount of excitations with 
energies significantly below our deduced Ag. However, 
until high quality single crystals become available it is 
impossible to rule out extrinsic origins of this term (for 
instance from weak links at the surface of the grains). 
The fit improves if two gaps are used (for instance with 
Al ^ 2A2 and the lower gap having a weight of 5% of 
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FIG. 2. Temperature dependence of the superfluid density 
p{T) = \{Qf/\{Tf for sample C (solid line). The dashed 
lines are the predictions of the anisotropic gap model (see 
text) 

the larger gap), but this leaves the value of the larger gap 
unchanged (within the error) from our original estimate. 

As mentioned above, several other authors have 
claimed that A(r) in MgB2 does not follow a simple expo- 
nential T dependence but rather a power law dependence 
with an exponent close to 2. In Fig. |^ we show a fit 
to data. Clearly, this fit is much worse than the BCS de- 
pendence [Eq. (^)]. We have also tried to fit other forms 
such as AT + BT"^ or AT^/(T + T*) but both fit the 
data much worse than Eq. (g) (the latter of these two 
equations is the form followed by a dirty d-wave super- 
conductor (l^). The key differences between our data 
and previous studies is our very high resolution, low base 
temperatures and very low excitation fields. 

As described above, for sample C we are able to de- 
duce the absolute values of A(T). In Fig. || we show 
the temperature dependence of the calculated superfluid 
density [p = A(0)/A(r)]^. Our data give a value of 
A(0) — 160 ± 20nm which is in good agreement with 
other estimates although is somewhat longer than 

that deduced from /iSR data 0. The uncertainty in A(0) 
does not make any significant difference to the tempera- 
ture dependence of p{T). The fit to the low temperature 
behavior of p{T) (Fig. |) gives values of Aq in the range 
28 K to 33 K as T^l^ was varied from 10 K to 15 K, 
in good agreement with those extracted from the fits to 
Ax(T) in Fig. ^ From this data we are also able to cal- 
culate the constant Ag in Eq. (^. We find that, unlike 
Aq, Aq depends quite strongly on the fitting range - in- 
creasing from 57 nm to 90 nm as T^ax varies from 
10 K to 15 K. This is however, significantly below the 
BCS value of A(0)= 160 nm. 

MgB2 has a anisotropic structure, with Mg atoms 
sandwiched between planar, hexagonal boron rings. It 
is therefore expected that there will be some anisotropy 
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FIG. 3. Low temperature behavior of p{T) for sample C 
along with a BCS fit (Eq. (||)). 

between the a5-plane and the c-axis responses. As our 
measurements have been performed on unaligned poly- 
crystalline samples our measurements of A(r) will be a 
mixture of the a6-plane and c-axis penetration depths (Ac 
and Xab)- Recent results on aligned crystallites and sin- 
gle crystals have shown a significant anisotropy in Hc2 
which implies an anisotropy in the coherence lengths, 
7 = S,ab/S,c- In an anisotropic Ginzburg-Landau the- 
ory this implies a similar anisotropy in A, 7 = Xc/Xab- 
There is some disagreement about the magnitude of this 
anistropy, Refs. [^,^ give 7 — 1.7 ± 0.1 and Refs. 
p2| , p3t give 7 = 2.6 ± 0.1. As far as we know the is 
no general solution for the moment of a sphere when A is 
anisotropic, however solutions do exist in the limits A 3> r 
Q and A ^ r Fortunately, these two limits give 

similar results for 7 ;< 3. To within ±10% the effective 
A(0) equals 1.2 Xab or 1.5 Xab for 7=1.6 and 2.6 respec- 
tively. Our value of A(0) then implies that AQf,=130nm 
and Ac=210nm or Aofc — 110 nm and Ac — 280 nm, for 
the two values of 7 respectively. 

Although the overall temperature dependence of A de- 
duced from our measurements is consistent with the a s- 
wave BCS picture, the value of the gap obtained is some- 
what smaller than the usual BCS weak coupling value 
(A/ksTc = 1.76). As pointed out by Haas and Maki 
||26| a simple explanation for this may be that the s-wave 
gap is anisotropic. In most conventional superconduc- 
tors the anisotropy in the pairing interaction is weak (~ 
10%). However, it has been proposed that the unique 
structure of MgB2 may produce a significantly enhanced 
anisotropy [p7| . 

Following Ref. |^ we model the anisotropy in A by 
A(z) = Ar(2:) where z — cos 6* {0 is the polar angle - 
with respect to the c direction), r(z) = (1 -t- az^)/(l -I- a) 
and a is the parameter which controls the anisotropy. To 
determine AA(T) in this model it is necessary to first 
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calculate A(T) self consistently from the gap equation 



1 



1 phu 



r{z)^dz 



JO (£:2 + A(z)2) 



and then AAaf,(T) from 



dEdz (3) 



A^(0) JoJa{z) dE (£:2_A(z)2)^ 

and similarly for \c{T) ( / is the fermi function). The 
temperature dependence of AA(r) measured on un- 
aligned polycrystalline samples can then be calculated 
by a suitable average of AAaf, and A Ac. For moderate 
to high anisotropy the low temperature behavior is dom- 
inated by AAaf, as in this model the gap is smallest in 
the plane. The resulting T dependence of A may then 
be fitted by Eq. (||) with an effective gap which is ap- 
proximately equal to the minimum gap in the model (the 
in-plane gap) . We find that in order to explain the exper- 
imentally determined value of A/Tc = 0.86 ±0.1 we need 
a — 2.2±0.4, or a total gap anisotropy (1 + a) = 3.2±0.4. 
Whether or not this very strong anisotropy may be ex- 
plained by a phonon interaction alone needs to be ad- 
dressed by microscopic theories. This anisotropy implies 
a significant difference in the temperature dependence 
of both A(r) and p{T) between the in-plane and out-of- 
plane directions. The different T dependencies of Pc{T) 
and pab{T) calculated from this model for a = 3 are 
shown in Fig. It should be noted although the detailed 
behavior will be different if strong coupling corrections 
and fermi velocity anisotropy are included, the qualita- 
tive difference between Pc(T) and Pab{T) will remain. As 
our measurements were made on unaligned polycrystals, 
our data will be an average of Pc{T) and Pab{T) (the 
weights depending on anisotropy in A(0)). Only a quali- 
tative comparison is possible, but it can be seen that the 
measured p{T) falls between Pc{T) and Pab{T). 

A different approach has been taken by Manske et al. 
psf , who have calculated the tunneling density of states 
from first principles within a standard Eliashberg the- 
ory. Their calculations suggest that the small value of 
Amm/^c is due to a pronounced low frequency depen- 
dence of A{uj) originating from a low-frequency phonon 
mode. The size of the minimum gap they calculate is in 
broad agreement with tunneling and our data. 

In conclusion we have presented a study of the pene- 
tration depth of polycrystalline MgB2 samples. In agree- 
ment with other direct probes of the symmetry of the 
order parameter we find that A(r) is well described by 
the usual s-wave BCS behavior, but with a gap that is 
significantly below the weak coupling BCS value. If this 
is to be ascribed to a simple anisotropy of the gap then 



we find that a fairly large anisotropy of around 3 is re- 
quired. Alternatively the recent calculation of Manske et 
al. j28| may also explain our data without the need for 
significant anisotropy in A(fc). Further studies on sin- 
gle crystal samples or aligned polycrystallites will help 
to clarify this question. 

We thank R.W. Giannetta, R. Prozorov, K. Maki, N.E. 
Hussey, J.F. Annett and J.R. Cooper for useful discus- 
sions, and P. Timms for providing the type B MgB2 sam- 
ples. 
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